package scu.maqiang.elasticity;

import scu.maqiang.fes.BVPType;
import scu.maqiang.mesh.Mesh2T3;
import scu.maqiang.mesh.T3Type;
import scu.maqiang.numeric.Direct;

public class PlainStress2DT31 {

	public static void main(String[] args) {

		Mesh2T3 mesh = new Mesh2T3().square2D(5, 5, T3Type.Right);
		
		Elasticity2DT32 elastic = new Elasticity2DT32(mesh);
		elastic.elastoStaticProblem(BVPType.PLANESTRESS, 
				new double[] {1.0e6, 0.25}, null, 
				new int[][] {new int[] {4}}, new Direct[] {Direct.All}, new double[][] {new double[] {0.0, 0.0}} , 
				new int[][] {new int[] {2}}, new double[][] {new double[] {-100.0, 0.0}}, 
				null, null, null, "elasticPlaneStress.dat");
	}

}
